Statistics of bubble rearrangement dynamics in a coarsening foam 



OO 
O 

o 

(N 



A.S. Gittings 1 and D.J. Durian 1 ' 2 

1 Department of Physics & Astronomy, University of California, Los Angeles, CA 90095, USA 

2 Department of Physics & Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA 

(Dated: July 15, 2008) 

We use speckle- visibility spectroscopy to measure the time-dependence of bubble rearrangement 
events that are driven by coarsening in an aqueous foam. This technique gives the time-trace for 
the average scattering site speed within a prescribed volume of the sample. Results are analyzed 
in terms of distributions of event times, event speeds, and event displacements. The distribution of 
rest times between successive events is also measured; comparison with diffusing-wave spectroscopy 
results shows that the spatial structure of a typical event consists of a core of only a few bubbles 
which undergo topology change plus a surrounding halo of bubbles which shift by an amount that 
decays to one wavelength at four to five bubbles away. No correlations are found between the 
durations, speeds, and rest times between successive events. 
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I. INTRODUCTION 

Aqueous foams consist of gas bubbles dispersed in a 
surfactant solution [H, H @]- With proper formulation, 
the soap films between neighboring bubbles are com- 
pletely stable against rupture. Therefore, the bubbles 
may be considered as packing units of fixed but de- 
formable volume. If the foam is sufficiently dry, such 
that the bubbles occupy a volume fraction greater than 
random-close packing for spheres, about 0.64, then the 
foam will exhibit a solid-like character with elastic mod- 
uli set by the ratio of surface tension to bubble size 0, Q • 
Since the bubbles are macroscopic, the corresponding en- 
ergy scale of modulus times bubble volume is many or- 
ders of magnitude larger than k^,T . Microscopically, the 
bubbles are therefore jammed, i.e. locked into a mechan- 
ically stable packing configuration where neighbors can- 
not change due to thermal motion [6|]. For bubbles to 
change neighbors, and for foams thereby to unjam and 
display a liquid-like character, requires an input of en- 
ergy. This may be supplied through an externally im- 
posed shear deformation. The energy for rearrangement 
may also be supplied through time evolution by coarsen- 
ing, in which gas slowly diffuses from smaller to larger 
bubbles so that the total interfacial area decreases Q- 
For very dry foams, rearrangements may be considered 
in terms of topology changes "of the first kind" , or Tl 
events, in the connected network of soap films PJ. 

The nature of microscopic bubble rearrangements and 
their connection to macroscopic foam rheology have thus 
been topics of fundamental interest. For bulk foams, 
Diffusing- Wave Spectrosco py (DWS) has emerged as an 
important tool 0, M (M [H 111 E H ■ This is a dynamic 
light scattering technique for opaque samples [la . Il6j . 
whereby far-field intensity fluctuations at an area com- 
parable to speckle size are measured and related to the 
absolute normalized electric field autocorrelation func- 
tion, gi (r) . The relative motion of scattering sites broad- 
ens the power-spectrum of the light, given by the Fourier 
transform of g\{r) exp(iw r) where co is the incident fre- 



quency, and causes the speckle pattern to fluctuate. Thus 
the functional form and decay rate of g\ (r) vs r give in- 
formation on the nature and the rate of the scattering site 
motion, respectively. For coarsening foams [3, Q, and for 
slowly sheared foams where the bubbles lock into place 
between rearrangements, [l3j . the contribution to g\(r) 
per scattering site is exponential in r with a decay time 
scale t set by the average time between successive rear- 
rangement events at a site. For rapidly sheared foams, 
where dissipative forces dominate surface tension forces 
so that the bubbles continuously flow and never lock into 
locally stable packing configurations, the contribution to 
g\ (r) per scattering site is exponential in r 2 with a decay 
time scale set by the strain rate and how long it takes for 
adjacent scattering sites to convect apart by the wave- 
length of light H3. 

While DWS thus reveals microscopic information 
through the shape and decay rate of the field autocorrela- 
tion, it cannot capture several key features of the bubble 
rearrangement dynamics. For example, it gives the av- 
erage time between events, but it cannot give the distri- 
bution. It cannot give the number of bubbles involved in 
an event. It cannot probe the bubble-size changes during 
coarsening and the resulting readjustment of the packing. 
Furthermore it cannot give the speed of bubble motion 
or the duration of the events. The latter is particularly 
important because it sets the crossover strainrate below 
which the foam exhibits discrete rearrangements and a 
solid-like rheology and above which it exhibits continuous 
microscopic flow and a liquid- like rheology [1J, [l7|, [3 • 
These limitations of DWS arise because the signal is ac- 
quired by time-averaging the speckle fluctuation statis- 
tics and also because the scattering volume must be large 
compared to the event volume. 

In this paper we apply a time-resolved dynamic 
light scattering technique called Speckle- Visibility Spec- 
troscopy (SVS) [Jjj|2C]] to a coarsening foam. The advan- 
tages over DWS are that the scattering volume need not 
be large compared to the event volume, and that speckle 
fluctuations are probed instantaneously by a single expo- 



sure of a CCD camera. From a time-series of such expo- 
sures, SVS allows the motion of individual rearrangement 
events to be examined directly. After reviewing details of 
this method, we present results on the detailed statistics 
of bubble rearrangements. This includes the full distri- 
bution of the time between successive events within the 
scattering volume, and by comparison with DWS results, 
an inference on average event size. It also includes the full 
distributions of event durations and peak bubble speeds. 



II. EXPERIMENTAL METHODS 

1 . Foam 

The experimental foaming system is Gillette Foamy 
Regular shaving cream, as in many prior studies dat- 
ing back to Refs. [a [9| and continuing as recently as 
Refs. [18|, |21(. The surfactants in aqueous solution are 
stearic acid and triethanolamine, plus other additives, 
giving a viscosity of fi — 1.8 cP and a surface tension of 
7 = 24 erg/cm 3 [1J|. Though no doubt the exact formu- 
lation evolves according to market influences, we notice 
no difference with any precedents [22J. Since this foam 
is produced by the aerosol method, the bubble sizes are 
smaller than possible with a mechanical mixing method. 
The advantage is that, therefore, it is possible to let the 
system age in-situ and coarsen significantly prior to data 
collection and thus to achieve a reproducible history- 
independent distribution of bubble sizes. Here we in- 
ject the foam into a rectangular glass cell with inner di- 
mensions 0.7 x 3 x 10 cm 3 , and seal. Data are collected 
throughout 4.4 — 7.1 hours after injection. During this 
period there is no drainage or film rupture, the liquid 
fraction is constant at about 0.08, and the average bubble 
diameter grows from 100 to 120 microns (roughly linear 
in time, since the interval is small, even though the ac- 
tual growth is a power-law with exponent near 1/2) [9(. 
The average time between coarsening-induced rearrange- 
ment events at a given scattering site, r , as measured by 
DWS P, grows from 40 to 55 seconds across this same 
age range. The average duration of events, r^, estimated 
by video of surface bubbles [a, EH 13 1 is °f order 0.1 sec- 
onds. Optically, bulk samples appear white and hence 
strongly scatter visible light with negligible absorption. 
The transport mean- free path, I*, equal to the typical 
step size in the random walk taken by photons as they 
diffuse within the sample, is 3.5 times the average bub- 
ble diameter [8|, (23|- The photon paths are not truly 
random walks, however, because photons are channeled 
with slight preference along the Plateau borders [2J| . 



2. Speckle-Visibility Spectroscopy 

The time-resolved dynamics of bubble rearrangements 
in the sample are probed by Speckle- Visibility Spec- 
troscopy (SVS). This method, reviewed in detail in 
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FIG. 1: (Color online) Example rearrangement events in a 
coarsening foam, illustrated by time traces of (a) raw pixel 
grayscale levels, (b) variance ratios computed from video data, 
and (c) bubble speeds computed from variance ratios. The 
foam age is to = 332 minutes, and the camera exposure time 
is T = 10 ms. The inset of (c) shows the bubble speed dis- 
tribution, averaged over age range 326 — 334 minutes, with 
the peak fit to a Gaussian. The threshold speed for discrim- 
inating events from noise is computed from the mode and 
standard deviation of the Gaussian fit as uth = Vmode + 3<r. 
This threshold is indicated by vertical and horizontal dotted 
lines in the inset and main plot of (c), respectively. The dura- 
tion of the events is denoted by Td, and the rest time between 
successive rearrangements is denoted by r r , as shown. 



Ref. [20] , has been applied to resolve the time-dependent 
motion of grains subject to periodic vibration [191 ] . of col- 
loids after cessation of shear ^25j , and of grains avalanch- 
ing down a heap [26|]. Closely related methods include 
laser-speckle ph otog rap hy J27l| and Time-Resolved Cor- 
relation (TRC) [23, [23, [3(|. Here we use the so-called 
"point-in / point-out" illumination / detection geometry, 
in which the volume of the sample probed by the detected 
light is not extremely large compared to the rearrange- 
ment event volume. This gives non-Gaussian field statis- 
tics for the fluctuations of a single speckle, which invali- 
dates DWS [3l[ but is ideal for SVS. Specifically, coherent 
light from a Nd:YAG laser, wavelength A = 532 nm, is fo- 
cused across a 3 mm diameter aperture formed by black 
tape attached and centered on the face of the sample 
cell. Photons which re-emerge from this same aperture 
are collected with a linescan CCD camera (1024 pixels, 
8 bits deep), placed such that the speckle size is compa- 
rable to the pixel size. There is no lens, only a line filter 



to eliminate ambient room light. The exposure time is 
set to T = 0.01 s, and the laser power is adjusted so that 
the average gray scale level is about 50. Video data are 
then collected in 8.5 minute segments, centered at ages 
{270, 300, 330, 360, 390, 420} minutes, and streamed 
to disk for post-processing. This is repeated four times 
starting with fresh foam, and the results are averaged to- 
gether for improved statistics; no differences are noticed 
between the individual runs. 

An example space-time plot for a 100 pixel subset of 
the video data is shown in Fig. [T{a) for a two second 
time window at age 332 minutes after injection. Static 
speckle is evident at the beginning and end of this win- 
dow as horizontal randomly-alternating bright and dark 
streaks; over these times the bubbles are nearly fixed in 
place. Two separate bubble rearrangement events hap- 
pen in between, centered near 0.7 and 1.1 seconds re- 
spectively, as evident by a scrambling of the speckle pat- 
tern. For proper choice of exposure time T, the speckle is 
easily visible between events but is smeared and washed 
out during events. In SVS motion is thus quantified by 
the variance of intensity levels, V^(T) ex (I 2 )t — (I) 2 , 
where (■ • -)t denotes the average over pixels exposed for 
duration T . Note that the average intensity is indepen- 
dent of T, so no subscript is placed on (I). The pro- 
portionality constant of VziT) is set by the laser inten- 
sity and the ratio of speckle to pixel size. Both effects 
may be neatly canceled by considering the variance ra- 
tio V2(m,T)/V2(T), where the numerator is found from 
a "synthetic exposure" equal to the sum of m successive 
images. As demonstrated in Fig. [TJb) for m = {2,4,8}, 
this variance ratio equals almost one when the speckle 
is static and decreases noticeably during rearrangement 
events - more so for larger m. 

The theory of SVS [IjJ, [2CJ may now be invoked to re- 
late the measured variance ratio to the underlying bub- 
ble motion. Specifically, if there is random ballistic speed 
between adjacent scattering sites with average speed v, 
then the instantaneous power spectrum is Lorentzian 
with linewidth 
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the corresponding normalized field autocorrelation is 
<7i(t) = cxp(— Tt), and the variance ratio is 
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where x = TT. From measurements of the variance ratio 
vs time, we invert this equation to deduce the linewidth 
and the average bubble speed v vs time. An example is 
shown in Fig.QJc). Note that the results for all three val- 
ues of m, give the same bubble speeds, validating our anal- 
ysis in terms of instantaneous linewidth T. However we 
emphasize one caveat: v represents the relative speed of 
scattering sites spatially averaged over the volume probed 
by the collected light. While diffusing photons may wan- 
der far into the sample, the primary sampling volume 
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FIG. 2: Autocorrelation of the velocity time-trace for a coars- 
ening foam over the age range 326 — 334 minutes. The ex- 
pected short and long time values based on the event and 
noise statistics of Fig. [3] are show by dotted lines, as labeled. 
The expected intermediate time value, labeled 1000 (nm/s) 2 , 
is given by /i V) 2 + /i(l - fi){vi){v ) + (1 - fi) 2 (v 2 }. 



for backscattered photons is approximately the aperture 
area times 2.5/* [32|. If a bubble rearrangement event is 
smaller than this scattering volume, then the speed de- 
duced from the linewidth will be an underestimate. We 
will present all data in terms of the nominal speeds given 
by Eq. ([T]). As will be discussed, this is accurate for the 
"noise" between events but is small by a factor of nearly 
10 4 for the events. 



III. TELEGRAPH ANALYSIS 

In this section we begin analysis of the SVS data on 
the time-dependence of the average bubble speed, v(t), in 
terms rearrangement event statistics. The first step is to 
identify the events, their durations Td, and the rest times 
T r between events. These quantities are independent of 
the connection in Eq. ([T]) between linewidth and bubble 
speed. Since v(t) is small but nonzero between events, as 
seen for example in Fig. [TJc), a threshold must be chosen 
for discriminating events from noise. For this it is helpful 
to consider the histogram of bubble speeds, as shown in 
the inset of Fig. [ljc) for the age range 326 — 334 min- 
utes. Usually the bubbles are not rearranging, so the 
speed distribution has a peak at small speeds represent- 
ing noise, and a long tail at high speeds representing 
rearrangement events. A reasonable prescription is to 



take the threshold as ^th 



3er where v mo de is 



the location of the peak in the speed distribution and a 
is the width of the peak given by fit to a half-Gaussian 
for v < Vmode- The result is indicated by a light-blue 
dashed line labeled i>th in both the main plot and inset 
°f Fig.QJc). Whenever v(t) rises above the threshold, we 
define the event duration r^ as the full-width at halfway 
between f mo dc and f max - Compound events are sepa- 
rated out whenever the signal remains above threshold 
but there is a minimum of depth greater than 3a below 



successive peaks. Rest times r r are denned as the interval 
between successive peaks. With these definitions, illus- 
trated graphically in Fig.QJc), the integral over an event 
is well-approximated by j) max T,j. And the fraction of time 
spent rearranging is given by the average duration and 
rest times as j\ = (Td}/(r r ). 

The integrity of these procedures may be tested using 
the autocorrelation of v(t), shown for example in Fig. [2] 
for the same 326 — 334 minute age range as Fig. Q] Fea- 
tures are to be compared with expectations based on con- 
sidering v(i) as a telegraph signal that switches between 
rest=0 and event=l levels with statistics given by the 
values of {(to), (r r ), f± — {Td)/{r r )} plus the mean and 
mean-squared speeds {(fo), (vq 2 ), (vi), (vi 2 )} computed 
respectively from P(v) with v below/above threshold. 
First note that the speed autocorrelation exhibits a two- 
step decay, with the first located near the average event 
duration (jd) and with the second located near the aver- 
age rest time (r r ). At short times, the autocorrelation is 
very close to /i(vi 2 ) + (l — fi)(vo 2 )- At long times, the au- 
tocorrelation is very close to [/i(i>i) + (l — /i)(fo)] 2 . At in- 
termediate times, between the two decays, there is an in- 

flection near / 1 2 ( Ul } 2 + / 1 (l-/ 1 )( Ul }( Uo ) + (l-/i) 2 («o 2 ). 
For a true telegraph signal this plateau would have been 
flatter. But the clear signature of all these features in 
the speed autocorrelation shows that we have achieved 
satisfactory identification of events. 

Scatter plots (not shown) of quantities {r^^+i, T r] j+i, 
w m ax,i+i} pertaining to event i + 1 vs quantities {r^^, 
i~ r ,i, ^maxji} pertaining to the previous event i reveal no 
evidence of correlations. Thus the switching between rest 
and event states appears to be Markovian, proceeding 
randomly at a history-independent rate. 



IV. EVOLUTION OF AVERAGES 
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FIG. 3: (Color online) Characteristic time and speed scales 
vs foam age. In (a) the solid diamonds represent the rest 
time between successive rearrangement events in the whole 
scattering volume and the solid squares represent the event 
durations, both averaged over 30 minute runs; the fraction of 
time the bubbles in the whole scattering volume are in motion 
is shown by the dashed curve /i . The time r between succes- 
sive events at one scattering cite was previously measured by 
diffusing-wave spectroscopy (DWS) and is shown by the solid 
curve. In (b) statistics for event speeds are shown by solid 
squares for the average, open squares for the standard devia- 
tion, and open pluses for the average of the maxima. Statistics 
for noise speeds are shown by solid diamonds for the average, 
open diamonds for the standard deviation, solid triangles for 
the mode, and dotted curve for the threshold; the latter two 
are from Gaussian fits to the full speed distribution as shown 
in the inset of Fig. QJc). The rate of change of the average 
bubble diameter, due to coarsening, was previously measured 
and is shown by the solid curve labeled "coarsening" . 



The average statistical quantities extracted from tele- 
graph analysis of the SVS speed signals are displayed vs 
foam age in Fig. [31 The average rest times and event du- 
rations are shown in the top plot, along with their ratio 
/i and the DWS results for the time r between succes- 
sive events at the same scattering site. There is good 
separation of time scales, r 3> (r r ) 3> (772). The dura- 
tions are constant for all ages studies; the average value 
is ((Trf)) = 0.44 s in accord with video observation of sur- 
face bubbles. The rest times grow longer with age, like r , 
but not dramatically; the average value is ((r r )) = 2.9 s. 
Note that these results for event and rest times are in- 
dependent of any underestimation in bubble speeds due 
to inapplicability of Eq. (fTJ) from interpreting localized 
motion as an average over the scattering volume. 

The key speed statistics are shown in the bottom plot 
of Fig. [3] For the events, this includes the average and 
standard deviation during events, and the average peak 
speed; these quantities are all approximately indepen- 
dent of foam age. As expected since the speed during 
an event rises and falls, the standard deviation is greater 
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FIG. 4: (Color online) Ratio of average rearrangement event 
diameter to bubble diameter vs foam age, based on SVS data 
for the average time (r r ) between events in the scattering 
volume and DWS data for the time to between events at a 
given scattering site. The average event size is 8.6 bubble 
diameters across for all foam ages (horizontal line). Error bars 
represent the range of results taking the scattering volume as 
aperture area times (2.5 ± 1.5)1* ■ 



than the average. Thus, in accord with a telegraph anal- 
ysis, the best quantity to consider is the peak speed; the 
average over all ages is ((v max )) = 120 ns/s. For the 
"noise" speed statistics, the plotted values include the 
average, the standard deviation, the mode, as well as 
the threshold used to discriminate events. These values 
are all comparable, are roughly constant [331 ] . and are 
much smaller than the event speeds. We also observe in 
Fig. [3] that the noise speed is ten times larger than the 
rate of change of the average bubble size, d(Db u bbie)/di, 
caused by coarsening. Therefore we infer that the noise 
in the SVS signal is due to the nonafhne relative motion 
of bubbles and plateau borders, continuously adjusting 
their positions in response to gradual size changes in or- 
der to optimize packing for locally-minimal surface area. 
While such motion maintains each bubble at the bottom 
of an energy well, the wells can eventually become unsta- 
ble and produce a sudden rearrangement event. 

The average time (r r ) between events in the scatter- 
ing volume can be understood in terms of the number 
of bubbles in an event. Since r r involves the whole 
scattering volume and the DWS time r involves only 
a single event volume, the value of T r is smaller than 
r according to the ratio of event to scattering volumes: 
T r = (Kvont / Scattering )t . Taking the scattering volume 
as aperture area times 2.51* [321 ]. and taking the transport 
mean free path I* as 3.5 times the average bubble diam- 
eter [8|, [23J , we thus deduce the ratio of event to bubble 
diameter and plot the value vs age in Fig. |U The results 
are nearly constant at -D even t/-Dbubbie = 8.6 independent 
of foam age. Though clearly finite, unlike avalanches 
which propagate out to the system size, this is still large 
compared to the bubble size. Video of surface rearrange- 
ments suggest that only a small core of neighboring bub- 
ble actually undergo topology change. Hence there is a 
halo of non-affine motion surrounding this "point-like" 
perturbation so that the bubbles maintain a locally op- 
timal packing. The extent of this reaction decays with 
distance, and falls below one wavelength at about 8.6 
bubble diameters. Fluctuations in the TRC signal are 
also related to event size, though the precise connection 
and the resulting event size were not reported [3J]. See 
Ref. [35| for an explicit discussion of how to deduce the 
size of dynamic heterogeneities from the variance in the 
decay of an autocorrelation. 



V. TIME DISTRIBUTIONS 

One advantage of SVS and the telegraph analysis is 
that the full distribution of all times and speeds are 
easily tabulated. For example, in Fig. [5] we display 
the distribution P(r r ) of rest times between events on 
scmilogarithmic axes. Since the average values consid- 
ered above are nearly independent of time, ((r r )) = 2.9 s, 
we combine together the results for all ages. The results 
exhibit a form which is nearly exponential, as demon- 
strated by comparison with the straight line representing 
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FIG. 5: Probability distribution for the rest time between 
successive bubble rearrangements in the scattering volume, 
averaged over all ages. For all events combined, the average 
is indicated by {{r r )) — 2.9 s. The dashed line represents an 
exponential distribution with this average. 
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FIG. 6: Probability distribution for the duration of individual 
events, averaged over all ages. The dashed curve represents 
the best fit to the sum of two Gaussians, with respective peaks 
at 0.29 and 0.54 s, with widths of 0.11 and 0.28 s, and with 
weights of 0.55 and 0.45. For all events combined, the average 
is indicated by ((t,j)) = 0.44 s. 



exp[— T r /(2.9 s)]/(2.9 s). This provides further support 
for the Markovian nature of the rearrangements noted 
above. Rearrangements occur seemingly at random with 
a given rate, just like the radioactive decay of unstable 
nuclei. The new packing states just after rearrangement 
must have a range of well depths. 

Next the full distribution -P(id) of the duration of 
events is displayed in Fig. [6] Again, since the averages 
are nearly constant, ((r^)} = 0.44 s, we combine the 
results for all ages. The form of the distribution rises 
abruptly from zero, peaks narrowly between 0.2 s and 
0.4 s, then decays more gradually out to about 1.2 s. Be- 
cause the rest times are exponentially distributed with an 
average only about ten times longer than the event du- 
ration, and because the scattering volume is larger than 
the event volume by the factor t /(t t ) = 17, it cannot be 
uncommon for two unrelated events to occur simultane- 
ously. Therefore we fit P{Td) to a sum of two Gaussians. 
The result shown by a dashed curve describes the be- 
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FIG. 7: Probability distribution for the event displacemen 
WmaiTd, averaged over all foam ages. The dashed line rep: 
sents fit to a power-law tail with exponent as labeled. For 
events combined the average is indicated by {{vmaxTd}}- 
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FIG. 8: Probability distribution for peak event speeds, aver- 
aged over all ages. The dashed line represents fit to a power- 
law tail with exponent —5/2 as labeled. For all events com- 
bined the average is indicated by {(t> ma x)). 



havior fairly well. The first peak is centered at 0.29 s, 
has width 0.11 s, and weight 0.55; this would represent 
singly-occurring events. The second peak is centered at 
0.54 s, has width 0.28 s, and weight 0.45; this would rep- 
resent multiple events. 

The average duration of events is a crucial microscopic 
quantity for understanding foam rheology, because dur- 
ing flow the nature of bubble motion changes dramati- 
cally according to whether the shear rate is high or low 
compared to the reciprocal of event duration [1J, [l?], [l8| • 
We note that the value found here, ((t<j)} = 0.44 s, is 
comparable to the relaxation time of a single Plateau 
border observed in a quasi-2D foam made with sodium- 
dodecylsulfate [361]. This time scale was argued to set 
by the sum of surface shear and dilatational viscosities 
divided by surface tension, independent of bubble size 
and bulk shear viscosity. In support, note in Fig. [3] that 
the event durations are independent of age even though 
the bubble radii and the time between events both grow 
noticeably. 
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FIG. 9: Probability density map for the speed vs time during 
an event; the color code for most to least probable is yellow, 
orange, red, maroon, brown, black. The speed axis is offset by 
the mode of the speed distribution and normalized by peak 
speed for each event; the time axis is offset by the time at 
which the speed is maximum and normalized by the event 
duration of each event. Only events in the power-law tail of 
the peak speed distribution of Fig. [8] are included. 



VI. MOTION DISTRIBUTIONS 

In this section we consider the nature of the bubble mo- 
tion during events. By the definition of event durations 
as the full-width half-max of peaks in the v(t) signal, 
the displacement due to an event is well-approximated 
by fmaxTd- The distribution of all such displacements, 
for all runs and all ages, is plotted on logarithmic axes 
in Fig. [7] It exhibits a peak around 12 nm, an average 
value of ((wmaxTd)) = 32 nm, and a long power-law tail 
with exponent —3. It is difficult to precisely interpret 
these results because, as noted above, the connection be- 
tween the SVS linewidth and bubble speed assumes that 
each scattering site has the same dynamics whereas the 
actual rearrangements are spatially localized. Further- 
more, the bubble motion itself is heterogenous, consist- 
ing of a core undergoing neighbor-switching and a halo 
undergoing nonaffine adjustment at fixed topology that 
decays in amplitude away from the core. Both effects 
may contribute to causing the apparent power-law tail. 
And certainly both effects contribute to making v(t) an 
underestimate of the bubble speed in the core. If we sup- 
pose that v(t) is dominated by the core motion, then we 
can correct for the underestimate by multiplying by (1) 
the ratio t /{tcl) ~ 50/0.44 of scattering to event volume 
and (2) the ratio (8.6/2) 3 of halo to core volume; this 
product equals 9000, the number of bubbles in scattering 
volume. The peak in P(v max Td) would then shift by a 
factor of 9000 from 12 nm to 110 /im. This is compara- 
ble to the average bubble diameter, which is the expected 
scale of displacement in the core. 



Further support for the above correction to the in- 
ferred speeds comes from consideration of the distribu- 
tion of peak bubble speeds, shown on logarithmic axes 
Fig. [SI Again this distribution is a compilation of re- 
sults for all runs and all ages. It exhibits a peak around 
30 nm/s, an average value of ((v ma x)} = 120 nm/s, and 
a long power-law tail with exponent —5/2. Applying 
the same factor of 9000 shifts the peak in the maxi- 
mum speed distribution from 30 nm/s to 270 jtxm/s. This 
value is in good agreement with an estimate of speed 
as equal to bubble diameter divided by event duration, 
v = (-DbubbieV^d) = 250 /im/s. Thus the analysis is 
self-consistent. 

Lastly we inspect the typical profile of speed vs time 
during an event in Fig.[3J Since there are broad distribu- 
tions of event durations and peak speeds, we normalize 
each event separately according to these characteristics 
and compile the results. To reduce the scatter, we also 
subtract the mode of the speed distribution, and include 
only events in the power-law tail of the peak speed dis- 
tribution, which, according to Fig. [5J begins at nominal 
speed of 50 nm/s. The color coding of bins in Fig. [9] is 
such that the most probable event shape is traced out 
along the crest of the yellow region. This reveals that 
events begin and end rather gradually, but speed up and 
slow down dramatically on either side of the peak. The 
event shapes are nearly symmetrical around the peak, 
with the rise being only somewhat faster than the decay. 
The first event in Fig. Q] is typical in its sharpness, while 
the second is more typical in its slight asymmetry. It 
would be interesting to know where the actual topology 
change occurs in relation to the peak speed. Perhaps the 
bubble accelerate toward a neighbor change which hap- 
pens at peak speed, or perhaps the bubbles remain at 
rest until motion is initiated by a neighbor change. 



VII. CONCLUSION 

For a very large scattering volume of foam, for instance 
with a "plane-in / plane-out" illumination/detection ge- 
ometry, multiple rearrangement events are present and 
each causes only slight speckle change since most light 
paths are unaffected. Thus the speckle fluctuations are 
nearly continuous and exhibit Gaussian field statistics at 
a single far- field detection site [3lfl . This is the usual re- 
quirement for intensity-correlation spectroscopy versions 
of dynamic light scattering like Diffusing- Wave Spec- 
troscopy (DWS). Under these conditions analysis of DWS 
data yields the rate of rearrangement events for coarsen- 
ing and slowly sheared foams, or the bubble speed for 
rapidly sheared foams. If the scattering volume is small 
enough so that often no rearrangements are present, and 
so that many light paths are affected by single events, 
then non-Gaussian field statistics in time-averaging at a 
single detection site will prevent application of DWS [3l[ . 
Under these conditions, there is extra information that 
can be accessed either by higher-order time-averaged in- 



tensity correlations [3l|, |37f or by time-resolved multi- 
speckle methods that ensemble-average over the speckle 
pattern. 

In this paper we probed coarsening-induced bubble 
dynamics in a small scattering volume with Speckle- 
Visibility Spectroscopy (SVS), which is a particularly 
simple and fast multispeckle dynamic light scattering 
method. The visibility of the speckle pattern in a sin- 
gle exposure of a CCD camera was repeatedly measured, 
then converted to a time trace of the linewidth of the scat- 
tered light. This signal exhibited a sequence of spikes rep- 
resenting individual discrete rearrangement events. Thus 
it was straightforward to identify the duration of events 
and of the rest times between events, which we analyzed 
in terms of both moments and full distributions. From 
lack of history dependence and from the exponential dis- 
tribution of rest times, we concluded that rearrangement 
events are Markovian. From comparison of rest times 
with previous DWS data for the time between events a 
a single scattering site, we concluded that the event size 
is almost nine bubble in diameter. Events consist of a 
small core of bubbles that switch neighbors, surrounded 
by a halo that shift at fixed topology. 

It was less straightforward to analyze the magnitude of 
the linewidth in terms of bubble speeds. Between events, 
the only motion is due to slight bubble size changes 
and gradual adjustment of bubble positions and Plateau 
borders to the changing packing conditions. This is 
stochastic but homogeneous both temporally and spa- 
tially throughout the scattering volume, and hence was 
readily analyzed in terms of bubble speeds. However, re- 
arrangements events are spatiotemporally heterogeneous 
so that the usual prescription for conversion of linewidth 
to speed gives a severe underestimate. We argued that 
the correction factor is set both by the ratio of scattering 
to event volume and the ratio of event to core volume. 
We showed that this gives both the correct total displace- 
ment within the core, and the correct speed within the 
core. In future studies, it would thus be a fine proce- 
dure to apply a scale factor to the measured linewidth so 
that the resulting typical displacement equals the bubble 
diameter. 

While the contributions here advance the understand- 
ing of both coarsening and rearrangement dynamics, 
and how they may be probed by time-resolved dynamic 
multiple-light scattering, there are many open lines of 
research. For example, it would be interesting to sys- 
tematically vary the scattering volume, both as a further 
check on the procedures presented here and also perhaps 
as a way of spatially-resolving the motion in terms of a 
diffuse-light tomography. Relatedly it would be interest- 
ing to consider the effect of the location of an event on 
the SVS signal. Also for example, it would be interesting 
to vary the aqueous solution in attempt to control and 
understand rearrangement rates and durations. And fi- 
nally, it would be interesting to understand the non-affinc 
shifts in response to coarsening or slow shear that ulti- 
mately lead to rearrangements. 
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